Astronomy & Astrophysics manuscript no. letter 


©ESO 2013 


March 8, 2013 





Letter to the Editor 



GRB 090618: a Candidate of a Neutron Star Gravitational Collapse 
to a Black Hole Induced by a Type Ib/c Supernova 

L. Izzo, Jorge A. Rueda, and R. Ruffini 

Dipartimento di Fisica and ICRA, Sapienza Universita di Roma, Rle Aldo Moro 5, 1-00185 Rome, Italy 
ICRANet, Rzza della Repubblica 10, 1-65122 Pescara, Italy 

ABSTRACT 

A novel concept has been recently proposed for explaining the temporal coincidence of some Gamma Ray Bursts (GRBs) with an 
associated Supernova (SN) in terms of the gravitational collapse of a neutron star (NS) to a Black Hole (BH), induced by a type 
Ib/c SN explosion. We apply these considerations to the exceptional case of GRB 090618, for which there is evidence of a SN ~ 10 
days after the GRB occurrence. We calculate the accretion rate and total accreted mass onto a NS from a SN Ib/c originated from a 
companion evolved star. It is shown that the NS reaches in a few seconds the critical mass and undergoes gravitational collapse to a 
BH leading to the emission of a GRB. We find for the mass of the NS companion, M NS , and for the SN core progenitor, M core , the 
following mass ranges: 1.8 < M N s /M Q < 2.1 and 3 < M corc /M < 8. Finally, we discuss the complementarity of these considerations 
to alternative processes explaining long and short GRBs. 
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The temporal coincidence of some Gamma Ray Bursts 
(GRBs) and their associated Supernovae (SNe) represent one 
of the outstanding and hottest topic in the Relativistic Physics 
and Astrophysics. We recall that the optical emission of a SN 
reaches its maximum about 15 days after the SN explosion (see 
e.g. Arnett 1996). This implies that the onset of SNe and the 
occurrence of the GRB can be almost simultaneous (see e.g. 
Della Valle 201 1; Hjorth & Bloom 201 1, for a recent review on 
GRB-SN systems). For these reasons the concept of the induced 
gravitational collapse (IGC) was initially proposed (Ruffini et al. 
2001) and further developed in (Ruffini et al. 2008). Recently we 
have theoretically developed the basic equations describing this 
concept (Rueda & Ruffini 2012). In (Rueda & Ruffini 2012) the 
GRB-SN connection has been explained in terms of the final 
stage of the evolution of a close binary system composed of an 
evolved star (likely a C+O star in the case of SN Ic) with a neu- 
tron star (NS) companion. The explosion of the SN Ib/c in this 
system leads to an accretion process onto the NS companion. 
The NS reaches the critical mass value in a few seconds, under- 
going gravitational collapse to a Black Hole (BH). The process 
of gravitational collapse to a BH leads to the emission of the 
GRB. We have evaluated the accretion rate onto the NS, and we 
have given the explicit expression of the accreted mass as a func- 
tion of the nature of the components and the binary parameters. 

We turn now to an explicit quantitative description of the pro- 
cess of IGC and compare the results with the observational data 
of GRB-SN. We apply in this letter these considerations to the 
specific case of GRB 090618 and its associated SN. We recall 
that GRB 090618 is one of the nearest (z = 0.54) and at the same 
time more energetic GRBs {E, so — 3 x 10 53 erg), ever observed. 
In fact, no beaming appear to be present from the data. These 
two important properties of GRB 090618 have produced an ex- 
ceptional set of data in view of the very fortunate joint observa- 
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tions of the major satellites observing in X and y rays: Fermi, 
Swift-BAT and XRT, AGILE, Suzaku-WAM, Konus-WIND and 
Coronas PHOTON RT-2. This has allowed to study GRB 090618 
in a very wide electromagnetic energy range, from 0.3 keV to 40 
MeV. There have been as well observations from the Optical tele- 
scopes: Hubble Space Telescope, Shane Telescope at the Lick 
Observatory and Palomar 60-inch Telescope, which led to the 
clear identification in the afterglow light curve of GRB 090618 
of a late ~ 10 days optical bump associated to the SN emission 
(Cano et al. 2011). 

We have recently shown that GRB 090618 (Izzo et al. 2012) 
is composed by two sharply different emission episodes (see 
Fig. 1). From a time-resolved spectral analysis, it turns out that 
the first episode, which lasts ~ 32 s in the rest frame, is char- 
acterized by a black body emission that evolves due to a de- 
creasing temperature with time (see Fig. 17 in Izzo et al. (2012)). 
Associated to the decreasing black body temperature, the ra- 
dius of the emitter has been found to increase with time (see 
Fig. 2 and Fig. 18 in Izzo et al. (2012)). With the knowledge 
of the evolution of the radius of the black body emitter, we 
find that it expands at non-relativistic velocities (see Eq. (1), be- 
low). Consequently, the first episode can not be associated to a 
GRB. Being prior to the GRB and therefore to the BH formation, 
this first episode emission has been temporally called Proto-BH, 
npaiTog, from the ancient Greek that means anteriority in space 
and time. 

We here identify the Proto-BH of the first episode as the first 
stages of the SN expansion. The black body emitting surface in 
the first episode evolves during the first ~ 32 s, as observed in 
the rest frame, following a power-law behavior 

r SN = at" , v SN = n — = ncrt'^ 1 , (1) 

where cr = 8.048 x 10 8 cm s _n , n « 3/5 as shown in Fig. 2, and 
\'sn = dr^/dt is the corresponding early SN velocity of the SN, 
so ~ 4 x 10 8 cm s _1 at the beginning of the expansion. 
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Fig. 1. Fermi GRB (Nal 8-1000 keV) light curve of GRB 090618. 
We explicitly indicate the two episodes of GRB 090618, identi- 
fied by Izzo et al. (2012). The first episode lasting 50 s in the 
observed frame corresponds to ~ 32 s in the rest frame of this 
source. 




Fig. 2. Evolution of the radius of the black body emitter given 
by Eq. (1) in the first episode of GRB 090618 (Izzo et al. 2012), 
here assumed as the radius of the SN. The time is here measured 
in the rest frame, so divided by a factor 1 + z, where z = 0.54, 
with respect to the one of the first episode given in Fig. 1. 

When the accreted mass onto the NS triggers the gravita- 
tional collapse of the NS to a BH, the authentic GRB emission 
is observed in the subsequent episode at t - to > 50 s (observer 
frame). The characteristics of GRB 090618 are shown in Table 
3 of Izzo et al. (2012) and we refer to that reference for further 
details on the GRB light curve and spectrum simulation. We find 
that the optically thick e + e~ plasma expands reaching a bulk 
Lorentz factor T ~ 500 at a radius of r IT = 1.4 x 10 14 cm. At 
the transparency of the e + e~ plasma the proper GRB emission 
has been observed at a temperature of 29 keV and the corre- 
sponding baryon load, B ~ 2 x 10~ 3 , has been obtained (see 
Sec. 5.2 in Izzo et al. (2012)). The baryon load is defined as 
B = Mgc 2 /E e to f , where Mb is the engulfed baryon mass from 
the progenitor remnant and E'J is the total energy of the e + e- 
plasma, which is assumed in the fit as the total observed GRB 
isotropic energy. We have obtained an average density of the cir- 
cumburst medium (CBM) of 0.6 particles/cm 3 , and each spike 
of the light curve corresponds to an interaction of the expand- 
ing plasma with clouds of radius as 10 15 cm and densities « 2 
particles/cm 3 (see Fig. 10 in Izzo et al. (2012), for details). 

We turn now to the details of the accretion process of the SN 
material onto the NS. The NS of initial mass M^s accretes mass 
from the SN ejecta at a rate given by (see Rueda & Ruffini 2012, 
for details) 



M acc (t) = np ej (t) 



(2GM NS ) 2 



rel,ej 



Pej(t) 



3M e ,(f) 



(0 



(2) 



where r 3 N (f) given by Eq. (1), M e j(f) = 
available mass to be accreted by the NS as a function of time, 



M ej , - M acc (t) is the 




Fig. 3. Time fQ CC since the SN explosion when the accretion pro- 
cess onto the NS starts as a function of the initial mass of the NS 
Mms and for selected values of the initial ejected mass M e j o, in 
the case of GRB 090618. 
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with M e jo the mass ejected in the SN. v re i,ej = -Jv, 
the velocity of the ejecta relative to the NS, where vsn is the SN 
ejecta velocity given by Eq. (1) and v rb = y/G(M cole + Mns)/o 
is the orbital velocity of the NS. Here M core is the mass of the 
SN core progenitor and a the binary separation. Hereafter we 
assume a = 9 x 10 9 cm, a value larger than the maximum dis- 
tance traveled by the SN material during the total time interval 
of Episode 1, At - 32 s, Ar ~ 7 x 10 9 cm (see Fig. 2). 

In the case when the accreted mass onto the NS is much 
smaller than the initial mass of the ejecta, i.e. M flcc /M e j,o << 1, 
the total accreted mass can be obtained from the analytic expres- 
sion given by Eq. (8) of (Rueda & Ruffini 2012), which in the 
case of GRB 090618 leads to 



MacM 



-f 



M acc (t)dt * (2GM m y 



15M ej , f 2/ 



8« 3 cr 6 Vl + kt 4 ' 5 



,(3) 



where k = v 2 . 1 (ncr) 2 and ?, a . cc is the time at which the accretion 

orb 1 v ' 

process starts, namely the time at which the SN ejecta reaches 
the NS capture region, R cap = 2GM^s / v rei ej ' so ^ or f — ^ we 
have M acc {t) = 0. The accretion process leads to the gravitational 
collapse of the NS to a BH when it reaches the critical mass 
value. Here we adopt the critical mass M a i t - 2.67M Q computed 
recently by Belvedere et al. (2012). The Eq. (3) is more accurate 
for massive NSs since the amount of mass needed to reach the 
critical mass by accretion is much smaller than M e j q. In general, 
the total accreted mass must be computed from the numerical 
integration of Eq. (2), which we present below for the case of 
GRB 090618. 

The occurrence of a GRB-SN event in the accretion induced 
collapse scenario is subjected to some specific conditions of the 
binary progenitor system such as a short binary separation and 
orbital period. The orbital period in the present case is 



P = 



4-n 2 a 3 



yG(M COK +M NS ) 



= 9.1 



M core + M h 

Mr, 



-1/2 



(4) 



We denote as Af acc the total time interval since the beginning 
of the SN ejecta expansion all the way up to the instant where 
the NS reaches the critical mass. In Fig. 4 we have plotted Af acc 
as a function of the initial NS mass, and for different masses 
of the SN core progenitor mass. The mass of the SN ejecta has 
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Fig. 4. Time interval At acc of the accretion process onto the NS 
as a function of initial NS mass M^s and for selected values of 
the SN core progenitor mass M core . The horizontal dashed line 
is the duration of the first episode of GRB 090618, At = 32.5 s, 
which constraints the duration of the time needed by the NS to 
reach the critical mass. The crossing points between the dashed 
horizontal line and the solid curves give the NSs with M^s that 
reach the critical mass in the time At. 



been assumed to be M e j o = M core - M rem , where M [em is the 
mass of the central compact remnant (NS) left by the SN explo- 
sion. Here we have assumed M core = (3-8)M at the epoch of 
the SN explosion, and M rem = 1.3M Q , following some of the 
type Ic SN progenitors studied in (Nomoto & Hashimoto 1988; 
Nomoto et al. 1994; Iwamoto et al. 1994). 

We can see from Fig. 4 that, in the case of GRB 090618, the 
mass of the NS companion that collapses to a BH, should be 
in the range 1.8 < M^s /M < 2.1 corresponding to the SN Ic 
progenitors 3 < M core /M < 8. 

The massive NS companion of the evolved star is in line 
with the binary scenario proposed by Ruffini et al. (2008). It 
is appropriate to illustrate a more accurate evolutionary sce- 
nario that leads to the occurrence of a GRB-SN event follow- 
ing the works of Nomoto & Hashimoto (1988); Nomoto et al. 
(1994); Iwamoto et al. (1994). The binary system is expected 
to evolve as follows: the initial binary system is composed of 
main-sequence stars 1 and 2 with a mass ratio Mo /Mi > 0.4. 
The initial mass of the star 1 is expected to be M\ > 11M Q , 
evolving through a SN event leaving a NS. The star 2, now 
with M2 > 11M after some almost conservative mass transfer, 
evolves filling its Roche lobe. It then starts a spiral in of the NS 
into the envelope of the star 2. Consequently, it will be composed 
of a Helium star and a NS in close orbit. The Helium star ex- 
pands filling its Roche lobe and a non-conservative mass transfer 
to the NS, takes place. This scenario naturally leads to a binary 
system composed of a C+O star and a massive NS, which can 
originate the process of IGC that we have proposed (Ruffini et al. 
2008; Rueda & Ruffini 2012). 

These results are also in agreement with the well under- 
stood Ib/c nature of the SN associated with GRBs. The most 
likely explanation for SN Ib/c, which lack H and He in their 
spectra, is that the SN progenitor star is in a binary system 
with a NS, see also Nomoto & Hashimoto (1988); Nomoto et al. 
(1994); Iwamoto et al. (1994) and also Tutukov & Fedorova 
(2007); Chevalier (2012) 

In the current literature, the GRB-SN connection has been 
at time assumed to originate from a specially violent SN pro- 
cess called hypernova or a collapsar (see e.g. Woosley & Bloom 



2006, and references therein). Both of these possibilities imply 
a very dense and strong wind-like CBM structure. This dense 
medium is in contrast with the CBM density found in most 
GRBs (see e.g. Ruffini 2011). Specifically, for GRB 090618 the 
average CBM density, («), inferred from the analysis of the after- 
glow has been shown to be of ~ 1 particle/cm 3 (see e.g. Fig. 10 
and Table 3 in Izzoetal. (2012)). The Lorentz Gamma factor 
of the electron-positron plasma has been estimated to be T = 
495+40, corresponding to a baryon load = (1.98+0.15)x 10~ 3 , 
(see Table 3 in Izzo et al. 2012). Such an ultra-relativistic expan- 
sion can only occur if the amount of engulfed baryon matter, 
quantifiable through the baryon load parameter B does not ex- 
ceed the critical value B ~ 10~ 2 (Ruffini et al. 2000). In view 
of the spherical symmetry of GRB 090618 this limit would be 
largely violated in the collapsar model. 

It is also interesting to compare the results on the IGC of a 
NS to a BH by a type Ib/c supernova Rueda & Ruffini (2012) 
with the results of Chevalier (Chevalier 1989) on the accretion 
of a supernova material by the central NS generated by the su- 
pernova. A total accreted mass of up to 0.1 M in a time of a few 
hours was obtained there for normal type II SN. Thus, a simi- 
lar amount of mass can be accreted in the two cases but in the 
latter the accretion occurs in a longer time. For reaching a high 
accretion rate of the inner SN material onto the central NS, there 
is the necessity of a mechanism that helps to increase the den- 
sity of the NS surrounding layers, which is decreasing due to the 
expansion after being unbound by the SN explosion. Chevalier 
(1989) analyzed the possibility of having a reverse shock wave 
that makes that work while it moves back through the SN core. 
The reverse shock is formed in the interaction of the mantle gas 
with the low density envelope. The time scale of the accretion 
process is thus determined by the time the reverse shock takes to 
reach the vicinity of the central newly born NS, which is of a few 
hours in the case of SN II progenitors. However, the existence of 
a low density outer envelope, e.g. H and He outer layers, is es- 
sential for the strength of the reverse shock. Fall back accretion 
onto the central NS is expected to be relevant only in SN II but 
not in SN Ic as the ones associated to GRBs, where H and He 
are absent. 

The argument presented in this letter naturally explain the 
sequence of events: SN explosion - IGC- BH formation - GRB 
emission. Correspondingly, the accretion of the material ejected 
by the SN into the nearby NS of the IGC model presented here 
occurs almost instantaneously. In fact, for the SN expansion pa- 
rameters obtained from the observations of Episode 1 in GRB 
090618 (see Eq. (1) and Fig. 2), the accretion of the SN material 
onto the nearby NS occurs in a few seconds (see Figs. 3 and 4). 
The binary parameters are such that the ejecta density does not 
decrease too much (from 10 6 to ~ 10 4 g cirT 3 ) before reaching 
the capture region of the NS, leading to a high accretion rate. As 
pointed out by Chevalier (1989), radiative diffusion will lower 
the accretion rate up to the Eddington limit (and then to even 
lower rates) when the trapping radius of the radiation in the flow 
r„. = KM acc l{4nc) (see Eq. (3.24) in Chevalier 1989), where k 
is the opacity, is equal to the Bondi radius r# = GM^s /v 2 el e •, 
the gravitational capture radius. The radius r lr is located where 
the outward diffusion luminosity is equal to the inward convec- 
tive luminosity. It can be checked that for the parameters of our 
system given by Eqs. (l)-(3), the equality r, r = rg occurs in a 
characteristic time ~ 200 days, where we have used k = 0.2 cm 2 
g -1 . Thus, this regime is not reached in the present case since 
the NS is brought to its critical mass just in a few seconds. In 
the case analyzed by Chevalier (1989), it happens in a time ~ 8 
days. 
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It is also appropriate, before concluding, to consider possible 
astrophysical systems which might have been the precursor of 
the GRB 090618 and its associated SN, as well as their possible 
final outcome. 

As for the precursor, a natural candidate are binary X- 
ray sources formed by a NS in a binary system with 
a main sequence companion, where the masses of the 
NS are observed to be close to the canonical value 
1.4M , e.g. Cen X-3 and Her X-l (Schreier et al. 1972; 
Wilson 1972; Tananbaum et al. 1972; Leach & Ruffini 1973; 
Davidson & Ostriker 1973; Gursky & Ruffini 1975); see also 
Rawls et al. (2011), for up-to-date estimates of the mass of the 
NSs in the X-ray binaries Vela X-l, 4U 1538-52, SMC X-l, 
LMC X-4, Cen X-3, and Her X-l. 

As for the possible final outcome, let us turn to the final 
state after the explosion of the SN and the trigger of the GRB. 
There is still a possibility of forming a NS-BH binary system. 
The outcome of the SN explosion leads necessarily to the forma- 
tion of a NS: evidence can be acquired by the observation of a 
characteristic late time X-ray emission (called URCA sources, 
see Ruffini et al. (2005)). This emission has been interpreted 
as originated from the young (f ~ 1 minute-(10-100) years), 
hot (T ~ 10 7 -10 8 K) NS, which we have called neo-NS (see 
Negreiros et al. 2012, for details). This emission can indeed be 
observed in the afterglow light curves of GRB 090618 (see e.g. 
Izzo et al. 2012) and GRB 101023 (see e.g. Penacchioni et al. 
2012). If the neo-NS and the BH are gravitationally bound after 
the GRB-SN emission, this new binary system might lead itself 
to a new merging process to a single BH. In this case the sys- 
tem could originate a further GRB with a possibly predominant 
emission in gravitational waves. 

In conclusion, the IGC binary scenario applied here to the 
specific case of GRB 090618 naturally leads to the understand- 
ing of the energetics and the temporal coincidence of SN and 
GRBs, as well as their astrophysical scenario and makes new pre- 
dictions on the final outcomes. It originates from a binary system 
composed of an evolved core and a NS. It is clear, however, that 
such GRB and their associated SN form a special class of long 
GRBs and of SNe Ib/c. There are in fact SNe Ib/c not associated 
to a GRB, e.g. SN 19941 (Immleret al. 2002) and SN 2002ap 
(Soria et al. 2004). Their observations refer to late phases of the 
SN evolution typically after ~ 15-20 days after the original col- 
lapse process. The existing descriptions of these late phases after 
15-20 days from the original explosion make use of a Sedov type 
behavior r oc f 2 / 5 (Sedov 1946, 1959). In the present case of the 
IGC we are presenting here for the first time, the first ~ 30 s of 
the very early evolution of a SN Ib/c associated to a GRB (see 
Eq. (1) and Fig. 2). The energetic of this SN Ib/c, as shown from 
Episode 1, appears to be much higher than the ones of the usual 
SNe Ib/c not associated to GRBs, Ei SOt Epa K 10 52 erg (Izzo et al. 
2012). The reason of this marked difference is certainly due to 
the accretion process during a SN explosion into the compan- 
ion NS and consequent gravitational collapse of the NS to a BH. 
The description of this challenging process although clear from 
a global energetic point of view, has still to be theoretically ex- 
plored in details and certainly does not present any relation to 
the Sedov type solution. 

From the point of view of the long GRBs, it is clear that the 
IGC concept leads to the formation of a BH of typical mass just 
larger than the critical mass of NS. 

There are of course GRBs originating in binary-NS systems 
and clearly do not have an associated SN, e.g. GRB 090227B 
(Muccino et al. 2012). These binary systems clearly depart from 
the one studied in (Rueda & Ruffini 2012) and in this Letter. For 



selected masses of the NS, such binary system leads to a short 
GRB with very low value of the baryon load, corresponding to 
genuine short GRBs. There are also systems as GRB 060614 
which do not show a SN (Delia Valle et al. 2006; Gal- Yam et al. 
2006) and are expected to be originating possibly from a yet 
different binary system formed by a white dwarf and a NS 
(Caito et al. 2009). 

Finally, it is appropriate to call on a selection effect in the 
study of IGC scenario. Only for systems with cosmological red- 
shift z < 1 the current optical instrumentation allows the obser- 
vation of the related SN Ib/c. A particular challenging analysis is 
the one of the system GRB 101023 (Penacchioni et al. 2012) in 
which the SN is not detectable but the IGC nature of the source 
is clearly recognized by the two different episodes in the GRB 
sources and the spectral features of the first episode. 
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